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Abstract 

In this paper, we discuss the time-space Caputo-Riesz fractional diffusion equation with 
variable coefficients on a finite domain. The finite difference schemes for this equation are 
provided. We theoretically prove and numerically verify that the implicit finite difference 
scheme is unconditionally stable (the explicit scheme is conditionally stable with the 
stability condition + < C) and 2nd order convergent in space direction, and 

(2 — 7)-th order convergent in time direction, where 7 G (0, 1]. 

Keywords: Time-space Caputo-Riesz fractional diffusion equation; Numerical stability; 
Convergence. 

1. Introduction 

Nowadays, fractional calculus has become popular in both the science and engineering 
societies. There are several nonequivalent definitions of fractional derivatives js], 0, 10 1. 



so 



14 



The Caputo derivative is most often used for time fractional derivative, and the Riemann- 
Liouville derivative and Griinwald-Letnikov derivative, the two fractional derivatives being 
equivalent if the functions performed are regular enough, are most frequently used for 
space fractional derivative. Some important progress has been made for numerically 
ving this kind of fractional PDEs by finite difference methods, e.g., see ,3, [HQ 

Another space fractional derivative having a vast majority of applications is the sym- 
metric fractional derivative, namely the Riesz fractional derivative, e.g., see 0, Il5|. Jiang 
et al analytically discuss the time-space fractional advection-diffusion equation with Riesz 
fractional derivative as space fractional derivative js]. Based on the shifted Griinwald ap- 
proximation strategy and the method of lines, Yang, Liu, and Turner numerically study 
the Riesz space fractional PDEs with two different fractional orders 1 < a < 2 and 
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< (3 < 1 [13[ . The explicit finite difference scheme for fractional Fokker-Planck equation 
with Riesz fractional derivative is discussed in ^ . With the desire of obtaining 2nd order 
convergence in the space discretization, here we further discuss the finite difference ap- 
proximations for two-dimensional time-space Caputo-Riesz fractional diffusion equation 
with variable coefficients in a finite domain, namely. 



u{x, y, 0) = uo{x, y) for (x, y) E 
^ u{x,y,t) = B{x,y,t) for {x,y,t) e dn x [0,T], 

in the domain Q = {xl,xr) x {yL^yn): < t < T, with the orders of the Riesz fractional 
derivative 1 < < 2 and the order of the Caputo fractional operator < 7 < 1; the 
function /(x, y, t) is a source term; and the variable coefficients c(x, y, t) > 0, d{x, y, t) > 0. 
The Riesz fractional derivative for nEN,n— l<u<n, in a. finite interval xl x < xr 



is defined as 



io| 



d''u{x,y,t) 



d\x\ 



where the coefficient /t^ 



1 



2cos{uTr/2) ' 



and 



1.2) 



XL 



D^^uix.y.t) 



T{n — u) dx'^ 



Qn rx 



n— i/— 1 



ui^,y,t)d^, 



;i.3) 



.D" u{x,y,t) 



xr 



X] 



:i.4) 



T(n — u) dx^ 

are the left and right Riemann-Liouville space fractional derivatives, respectively. The 
Caputo fractional derivative of order 7 G (0, 1] is defined by [8|, l9| 



^D;'uix,y,t) 



r(i-7).o 

du{x,y,t) 



du{x,y,r]) 



{t - 77) ^c??7, < 7 < 1, 



;i.5) 



dt 
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1. 



For the 2nd order discretization of the Riemann-Liouville fractional derivatives (11. Sp 
and (II. 4p . it has been detailedly discussed in [l| being the sequel of 11 1. This paper applies 
the discretization scheme to Riesz fractional derivative. The implicit and explicit finite 
difference schemes are designed. We theoretically prove that the implicit finite difference 
scheme is unconditionally stable and the stability condition of the explicit scheme is 
+ < C, confirming the general conclusions on the stability conditions of the 
explicit schemes for fractional PDEs The desired 2nd order convergence in space and 
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(2 — 7)-th order convergence in time of both implicit and explicit schemes are theoretically 
proved and numerically verified. 

The outline of this paper is as follows. In Section 2, we introduce the approximation 
of the Caputo fractional derivative and the 2nd order finite difference discretizations for 
the Riesz fractional derivatives, and derive the full discretization schemes of (11. ip . In 
Sections 3 and 4, the stability and convergence of the provided implicit and explicit finite 
difference schemes are analyzed, respectively. To show the effectiveness of the schemes, we 
perform the numerical experiments to verify the theoretical results in Section 5. Finally, 
we conclude the paper with some remarks in the last section. 

2. Derivation of the finite difference scheme 

We use two subsections to derive the full discretization schemes of (ll.ip . The first 
subsection introduces the approximation of the Caputo fractional derivative and the sec- 
ond order finite difference discretizations for the Riesz fractional derivatives in a finite 
domain. The second subsection gives the full discretization scheme (implicit scheme and 
explicit schemes) to the one-dimensional case of (II. ip and (II. ip itself, respectively. 

2.1. Discretizations for the Caputo and Riesz fractional derivatives 

Take the mesh points Xj = xl + iAx, i = 0,1, . . . , N^, yj = Vl + j^U, J = 0, 1, . . . , A^y 
and tk = kT,k = 0,1, ... , Nt, where Ax = {xr - xl)/N^, Ay = {yn - yL)/Ny, r = T/Nt, 
i.e.. Ax and Ay are the uniform space stepsizes in the corresponding directions, r the time 
stepsize. For v G (1,2), the left and right Riemann-Liouville space fractional derivatives 
(II. 3p and (ll.4p have the second order approximation operators S^^^^u^j and S^^^^u^j^j, 
respectively, given in a finite domain [ij], where u^^ denotes the approximated value of 
u{xi,yj,tk). 

The approximation operator of (II. 3p is defined by [jj 




(2.1) 



and there exists 



,D>(x,^/,t) = 5,,^,<^. + 0(Ax)^ 



(2.2) 



where 




m < i 



i - 1 



m = I, 



(2.3) 



= < 




m = i + 1, 

m > i + 1, 
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and 



(j-l)3--j2-(j-3 + z/), m = 0, 

-3,rn - >i (j - m + If''' - 2(j - mf-'' + {j -m- if-"", 1 < m < j - 1, 
1, m = j, 

with j = i — 1, i, i + 1. 

Analogously, the approximation operator of (11. 4p is described as 



S ■ ■- 



and it holds that 



with 



^ ' ^ ' m=i—l 



Q = < 



0, 

bi-i,i-i, 



m < i — 1, 
m = i — 1, 
m = i. 



and 



(2.4) 



(2.5) 



(2.6) 



1, m = j, 

bj,m = { (m - J + 1)3-" - 2(m - j)^-" + (m - j - 1)^-", j + 1 < m < - 1, 
i3-a-N^+ j){N^ - jf-" + (iV, - J - 1)3-", m = iV„ 
with j = i — + 1. 

Combining (12.21) and (12.51) . we obtain the approximation operator of the Riesz frac- 
tional derivative 



d''u{x,yj,tk) 



d\x\'' 



where 



V 

9i,m 



iu {5u,+x + + 0(A 



X 



r(4 - u)Ax 



m=0 



r(4 - z/)Ax^ 



m=0 



Vi ml 



m < i — 1 

Pi,i-i + = « - 1 

m = i 

m = i + 1 



a" 



m > 2 + 1. 



(2.7) 



(2.8) 
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Taking = 2, both Eq. (12. 2p and (12. 5 P reduce to the following form 

d'^u{xi, y, t) u{xi+i, y, t) - 2u{xi, y, t) + u{xi-i, y, t) 



+ C(Ax) 



The Caputo derivative in the time direction is discretized as j4| 



^ I > s=0 



"*"+^ du{xi, yj, 7]) dr] 

dr] {tk+i - r;)T 



k 



^ sr^ u{xi,yj,ts+i) -u{xi,yj,ts) f dr] 
r(i-7)t^ r k {tk+i-vV^ ^ 

k 

]- , ujxj, yj, tk+i-s) - u{xi, yj, 4-.) 2-7n 

I _ ^ T-7 + l'^ J) 



r(2-7)f:^ 



(2.9) 



where 7 G (0, 1), Is = {s + 1)^"^ - s^"^. 

When < 7 < 1, the time Caputo fractional derivative uses the information of the 
classical derivatives at all previous time levels (non-Markovian process). If 7 = 1, then 
Iq = 1, Is = 0, s > 0, it can be seen that by taking the limit 7 — 1 in (12. 9p . which gives 
the following equation 

du{x,y,tk+i) _ u{x,y,tk+i) - u{x,y,tk) 



dt 



0{t). 



Similarly, it is easy to get the one- dimensional case of (I2.ip -( !?]9|) . 



Remark 2.1 



Denoting f/" = [u^ p u^j, ■ ■ ■ , W^^-i j] , j = 0, 1, . . . , A^j^ and rewriting 
and ([22D as matrix forms 5„,+^.f/" = if/" + 61 and 5^^_Jj'' = BU'' + b2, respectively, 
then there exists A = B^. 

2.2. Implicit and explicit difference schemes 

Let us first consider the one-dimensional time-space Caputo-Riesz fractional diffusion 
equation 

^D7^(x,t) = c(x,t)^^|^ + /(x,t). (2.10) 
Using the one- dimensional case of (I2.ip - (l2.9p . we can write (I2.10p as 



7--7 ^ 



r(2 - 7) ^ 



^ ls[u{Xi, tk+l-s) - u{Xi, tk-s)] 



K^c[x,,t,^,)_ ^ g.^^^^^^t^^^^ + /(x.,4+0 + Oir'-' + (Ax)^), 



r(4-a)(Ax, 
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where i = 0,1,..., N^, k = 0,1,..., Nt. Assuming that = c{xi, 4), fj' = f{xi, 4), 
Ui^k+i = - ^r{A-a){Ax)'- /i = r(2 - 7)r^, we have 

k 



u{Xi, tk+i) = u{Xi, tk) - ^ ls[u{Xi, tk+l-s) - u{Xi, tk-s)] 

s=l 



(2.11) 



m=0 



where < Ct'{t^-^' + {Axf). 

Therefore, the imphcit difference scheme of fl2.10p has the following form 



J2 Is^'-' - u\-') + + (2.12) 

s=l m=0 

and it can be rewrote as 

(1 - ^i,l9t^u\ - UJi^i Y gtm^ln = U^i + f^fl, k = 0, 

N:, k-1 

(1 - u,,,+,gr,)u^+' - ^ = 5^(/, - /.+i)«r^ + hu^ + /x/f +\ > 0. 

m=0,m^i s=0 

(2.13) 

Taking 7 = 1, thus the implicit difference scheme fl2.13p reduces to the following 
equation 

(1 - a;,,+i^?,^,)«^i - J2 9t^<^' = u\ + /i/f +\ k > 0. 

Next, we examine the two-dimensional time-space Caputo-Riesz fractional diffusion 
equation (11. ip . According to (I2.ip - (l2.9p . then (II. ip can be recast as 



— — — ^/,[n(xi,?/j,tfc+i_,) -u{xi,yj,tk-s)] = - " 9?,mU{x^, Vj, 4+i) 

- y(4!"^^(l7/j i; + fi^^.yvt,^l) + 0{T^-- + (Ax)^ + (Ay)^), 

where i = 0, 1, . . . , A^^;, j = 0, 1, . . . , A^y, k = 0,1,. . .,Nt. Denoting c^j = c{xi, yj, tk), 
di,j = d{xi,yj,tk), ftj = f{xi,yj,tk) and 

r(2-7)r^^.c:-+\ _ r(2 - 7)r^/.^4+\ 



r(4-a)(Ax)" ' r(4-/3)(Ay)/5 ' 



6 



and /i = r(2 — 7)r''', we obtain 

k 
s=l 

/ (2-14) 

+ ^i,j,k+l 2^ ^*"mW(a;m, Vj, tk+l) + ^i,j,k+l 2^ 9j,m'^iXi, y^, tk+l) 
m=0 771=0 

+ fif{xi,yj,tk+i) + Ri;^\ 

where < Cr^(r2-^ + (Ax)^ + (Ay)^). 

Then we obtain the full discretization implicit difference scheme of fll.ll) as 

k Ny 
s=l m=0 m=0 

(2.15) 

and Eq. (I2.15P can be rewritten as 



(1 ~ ^'i,j,l9i,i - ^i,j,l9j,j)ulj - Uj'ij^i ^ 9?,m'^ln,j ~ X/ 9j,m'^],m 

~ ^i,j,k+l9?,i ~ '^i,j,k+l9jj)lJ'iJ^ ~ '^i,j,k+l 9i,m^^} ~ ^i,j,k+l 9 

m=0,m^i m=0,mj^j 

k-1 

= - /.+i)<7^ + + fift;\ k > 0. 

s=0 

(2.16) 

When 7 = 1, fl2.16p becomes 



(1 ~ ^i,j,k+l9i,i ~ ^i,j,k+l9j,j)Uij^ - ^i,j,k+l 9Lm'^'^,j ~ ^i,j,k+l 9 

Let u*^ = [ug, u^, . . . , u^J^, f*^ = [fo^ ff^, . . . , f^J^, where = [uf^g, Mf^^, ...,u'y 



and ff = [fl^Q, fl'i, . . . , fi^NyV- Then ( I2.16P can be written in the matrix form 



Avi- = u° + /if\ = 0, 

^•-1 (2.17) 



s=0 

where A is a (N^Ny) x (N^Ny) coefficient matrix. 
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Analogously, the explicit difference scheme of fl2.10p is 

No: 

u\ = (1 + a.^ogtM + ^^.0 9tmul + l^fl k = 0, 

m=0,m^i 
fc-1 

U-^' = (1 - /l + (T.^kQiM + J2ils - + + (T,,fc ^ + /i/f , /C> 0, 

s=l m=0,mj^i 

(2.18) 

where CTj^fc = ~ ^r{4-l){Ax)°' ^ ^'^'^ ^^-^^ explicit difference scheme of fll.ip can be 
expressed as 



(2.19) 



u\f = (1 - /i + + fc(7,^,)M^,. + J](/. - /s+i)mJ-^ + 

s=l 

^i,3yk / ^ 9i,m^m,j "I" '^i,j,fc / ^ 9j^ni^i,in ~l~ /^/ij' fc > U, 

, _ r(2-7)r^^ccfe^- „ _ r(2-7)r^K^dfc^ 

wnere CTj^y. r(4-a)(Ax)" ' ^jj.fc r(4-/3)(A?;)/3 • 

3. Stability analysis 

Now we perform the detailed stability analysis for the implicit and explicit schemes 
fl2.13p and (12.180 of the one dimensional case (I2.10p . and the implicit and explicit schemes 
(I2.16P and (I2.19P of the two dimensional case (ll.ip . First we introduce two lemmas on 
the properties of the coefficients of the discretized fractional operators. 
Lemma 3.1 

aa. 

Let 7 G (0, 1), then coefficients Is defined in ^2.9^) satisfy 



(1) >0,s = 0,l,...,A;. 

(2) \ = Iq > li ■ • ■ > 4, /fc as k ^ oo. 

(3) Cik'^ < (/fc)~^ < C2k'^, where Ci and C2 are constants. 

k 

(4) ^(/, - /,+i) + h+i = (1 - h) + Yi^s - h+i) + lk = l. 

s=0 s=l 

Lemma 3.2. The coefficients g^^, v G (1,2] defined in k2.8i) satisfy 
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(1) 0, gl^>Q{m^i)- 

(2) 5Z<„.<0 and -gl,>J2gl^. 

m.=0 m=0,mj^i 

Proof. According to Theorem B in the Appendix of and formulas (12. 3p and (12.61) . 
we obtain q^^ > for m > i + 1, p^^ > for m < i — 1, and p^j+i = <lii-i = 1, 
Pii = lii = — 4 + 2^""^. From the definition of g^,^ (12. 8p . it is easy to verify that (y'j'^j+i = 
g'j;.^^ = 7- 2^"" + i^^" > 0, ^J^i = -8 + 2"^-" < 0, and other gl^ are positive. And letting 

M be a constant in (12.71) . we can easily derive that Yl dim < 0' then —g^i > Y 9im- ^ 

m=0 m=0,mj^i 

Next, we use four subsections to strictly prove that both the implicit schemes (I2.13P 
and (I2.16P are unconditionally stable; the explicit scheme (I2.18P is stable under the condi- 



tion < C and the explicit scheme (I2.19P is stable under the condition + 



< 



C. 



3.1. Stability for one- dimensional implicit difference scheme 

Theorem 3.3. The implicit difference scheme ^2.13\l of the one- dimensional time-space 
Caputo-Riesz fractional diffusion equation Ii2.10\) with 0<7<l,l<«<2zs uncondi- 
tionally stable. 

Proof. Let (z = 0, 1, . . . , N^; k = 0,1, . . . , Nt) be the approximate solution of u^, which 
is the exact solution of the implicit scheme (I2.13p . Putting ef = — Ui, then from (I2.13p 
we obtain the following perturbation equation 

k-1 ^^-^^ 

(1 - )e^i - g^^^e'^' = ^(/, - U^)et' + ke^ k > 0. 

When 7 = 1, Eq. (13. ip can be written as 

(1 - i^.,k+igtM^' - ^^.'^•+1 E 3tm^'^' = 4. k > 0. 
Denoting = [en,e?, . . . , 1 and ||-E''^||oo = max lefl, then we use the mathematical 

0<i<Na: 

induction to prove the unconditional stability. For k = 0, supposing le^^^^l = \\E 
max le.^l, according to Lemma 3.2, we get 

0<i<N^ ' * ' 

\\E ||oo = l^igl ^ l^iol ~ ^to,l 9io,m\^io I ~ \^io \ ~ '^«0,l5'jo,io 1^*0 I ~ '^^O,! ^ ] 9io,m\^io 

m=0 m=0,m^io 



Ico 



a A 



< (1 - Wio,l<,io)l4l - ^^0,1 9Z,m\^m\ < |(1 - ^io,l9Z,iMo - ^»o,l 3io,m^ 

= le? |< 



|00' 



Assuming IIE^IL < ||^°||oo, A; = 1, 2, . . . , fc, and le'^+^l = I L = max |ef+^|, there 

" 0<i<N^ 

exists 

I loo - \% I S \% I -Wio.fc+l 2^fi'io,mFjo I 

m=0 

— l^io \ ~ ^io,k+igio,io\% \ - ^io,k+l 2^ 9io,m\^io I 

m=0,mj^io 

m=0,m^io 

< |(1 - Ui^^k+igta,ki)4o^ - ^io,k+l Y, 9io,m^^m^\ 
k-1 

=0 

k-1 

^k—s I 



k~s I 



I oo 



oo 



From Lemma 3.1, we obtain 

k-1 

WE^^'Wo. < (1 - ^l)ll^1loo + 4||i^°||oo + Y^ls - ls+l)\\E 

s=l 
k-1 

< (1 - /l)||^°||oo + /fc||^°||oo + Y^ls - 

s=l 

< ||-^°||oo- 

When 7 = 1, using similar idea, we can prove 

ll-E^'^^lloo < 1(1 - Wjo,fc+lfl'j",io)eio^^ - ^io,k+l X] 9io,m^^m'^ \ = l^io I - H"^"! 

3.2. Stability for one- dimensional explicit difference scheme 
Theorem 3.4. // 

r(4-a)(l -2-^) , ^ . ^ 

< -7-. — ^ < — ; 7-, ; — r, where Cmax = max ciXiAu), 



00 • 



□ 
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then the explicit difference scheme Ii2.18\) of the one- dimensional time-space Caputo-Riesz 
fractional diffusion equation Ii2.10\) with 0<7<1, 1 < a < 2 is stable. 

Proof. Under the above conditions, we obtain < —o'i j.g"'- < 2 — 2^^"' and l—li + cjikg"^ > 
0. Assuming that {i = 0,1, . . . , N^; k = 0,1, . . . , Nt) be the approximate solution of u^, 
which is the exact solution of the explicit scheme fl2.18p . Therefore the error = mJ^ — 
satisfies 



(3.2) 



e\ = (1 + (Tiflgli)e^i + ai,o 9trr/m^ k = 0, 

s=l rn=0,mj^i 

When 7 = 1, fl3.2l) becomes 

e^^ = (1 + a.,k9ti)e^ + a.,k Yl ^tm^'m, k > 0. 

Denoting E'' = [eg,ei, . . . and Hi^'^Hoo = '^^ use mathematical induction 

to prove the conditional stability. For k = 0, supposing \e}J = \\E^\\oo = q™-^^ 14 
according to Lemma 3.2, we get 



||-£^^||oo = ICiol = |(1 + Crio,05'j",io)ei'o + ^^^0,0 Y 9io,m(^m\ 

m=0,m^io 

,m I I 

m=0,m^io 

< (l + a.„,o<,J||i?°||oo + cr.o,o 5^ 9tJ\E'\ 
= \\E'\\o. + ^^o,oJ29ZJ\E'\\oo<\\E'^' 



oo 



oo • 



m=0 



Assuming ||-E'*'||oo < ||-E'°||oo, k = 1,2, . . . ,k, and denoting |ef+^| = ||£'''+^||oo = max |e 

0<i<Nx 

from Lemma 3.1, there exists 



fc+ii 



fc-l iV^ 

ll^'+'lloo = \ei^'\ = |(1 - /i + ^.o,;^<.o)4 + J2(^s - ls+i)ei-' + kel + a.,,u ^^c^d 

s=l m=0,mjtio 
fc-1 Afx 

< (1 - /i + cx^,,k9ZJ\el\ + Yi^s - Wi)|etri + ^f^Kl + ^*o,fc dZJ^tl 

s=l m=0,mj^io 
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k\ 



< (l-/0ll^1|oo + 5^(/.-Wl)||i^^'=-1|oo + 4||i?°||oo + a,,,,^<„J|E 

s=l m=0 

< (1 - /i)||i5;°||oo + J2(is - Wi)ll^°lloo + hWE"^' 



I oo 



s=l 



(1 - /l + ^(/, - Ui) + 4)||i?°||oo = 



If 7 = 1, using similar method, we can prove 

m=0,mj^io 

< (1 + (^io,kgZ,J\ei I + (Tio,k 3l ,m I I 

m=0,mj^io 

< (l + ^.o,<.o)ll^'lloo + ^.o,;^ E 

= \\E'\\^ + a,,^,J29lJ\E'\\oo < WE'W^ < \\E' 

m=0 



oo 



□ 



3.3. Stability for two-dimensional implicit difference scheme 

Theorem 3.5. The implicit difference scheme 1^2. 16\) of the two-dimensional time-space 
Caputo-Riesz fractional diffusion equation U.l\) with 0<7<1, 1 < a, f3 < 2 is uncondi- 
tionally stable. 

Proof. Let u^j {i = 0, 1, ... , N^; j = 0,1,..., A^^; k = 0,1,..., Nt) be the approximate 
solution of ufp which is the exact solution of the implicit scheme fl2.16p . Denoting that 
e^ j = u^ j — u'lj, from (12. 16^ we get the following perturbation equation 

(1 ^ ^i,j,k+l9i,i ~ ^i,j,k+l9jj)^ij ~ ^i,j,k+l 5^ 9i,m^m,j ~ ^i,j,k+l 9j,m^i,m ~ k = 0, 

Ny 

l-"- ^i,j,k+iyi,i ^i,j,k+iyj,jJ^i,j '^i,j,k+l yi,m^m,j ^ij,k+l 9j,m^i,m 
k-1 

= E(^s-WiK-' + /;^<,-, k>0. 

s=0 

(3.3) 
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When 7 = 1, (13.31) becomes 



No, Ny 
\^ ^i,j,k+iyi,i ^i,j,k+iyj,j)^i,j ^i,j,k+l yi,m^m,j ^i,j,k+l yj,m'^i,m ft. ^ U. 



Denote u'' = [Ug, u^', . . . , u^^]'^, u*' = [r^ 



1 , • • • , 



k IT 



and I IE* 



We prove the resuhs by mathematical induction. For = 0, supposing \(\ 



max le.^J 
I |E"^| Ifv^ = 



«0 JO I 



max 

0<j<Ar^,0<i<Ary 



there exists 



lE^I 



«0,J0 I 



m=0 



m=0 



"^O JO 



' fl" I _ \ ^ I 1 

iOjiOiiyiOiiol «OJol «0,jo,l / J yio,m\Ht 



«0 JO I 



Ny 



— " ''^ I ^ I — " ^ I 

'^joj'o,l5'jojbl^«Oj'ol '^ioj'o,! / J yjo,m\^k 

m=0,m^jo 

<(1-<J0,1^^ 



«0 JO I 



,/3 



«ojo,i^«o,«o '^'iojo, JO /I «0J0 

. I — r I ; 

mjo 



/3 1,1 I 



< 



^ "^IQ JO, 1^40,10 ioj'o.li'joj'O'' «( 



^«Oj'o,l / y fi'io.ml^ 
m=0,m^jo 

/3 



«o JO 



Ny 



«0J0,1 / ^ yio.m'^mjo «OJO,1 / ^ yjo^m'^M 



iO,m 



"^o JO 



< llE^llno. 



m=0,m^jo 



Supposing II E'^ I loo < ||E°||oo, k = 1, 2, . . . , /c, and |ef+/g, 
then 



IE 



fc+i| 



oo = max 

0<i<N^,0<j<Ny 



1 



IE 



fc+l| 



I / I fc+1 I _ / \^ n° Ip'^+M 

«oJol — I «oJol ^■tpj'o.fc+l / ^ ytn,m\'^i(),i()\ 

m=0 



Ny 



«OjO,fc+l / ^ yio.mFjoJol 
m=0 

^fc+l I 



fc+1 I 
io JO I 



«oJo,fc+iy«o,«ol iojol «oJo,fc+l / ^ y«o,mi-iojoi 

m=0,m^io 



\^k+l 



'^iojo.fc+l^iojol^^ojol '^joJo.fc+l ^io,m 

m=0,m^jo 



io Jo I 



— '^iojo,fc+lfi'jo,«o '^iojo,fc+lfi'iojo)l^«o'jol 

N^ 



Ny 



oo' 



«Oj'o,fc+l ^ y yio,m\-m,jo\ 
m=0,m^io 



ioj'o,fc+l / , ^]o,m\-io,in\ 
m=0,m^jo 
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< 



Ny 



fc+1 



m=0,m^io 



I 



s=0 

By Lemma 3.1, we get 



IE 



fc+ii 



< (l-^i)llEl 
<(l-/i)l|E°| 

< I I E*^ I Inn. 



m=0,m^jo 



fc-1 



I -I- / I I'F^l I 
I oo ~r '■fc 1 1 -•-^ 1 oo 



I / I IttOi I 

oo ~r '■fc -'-^ oo 



k—s I 



s=l 
fc-1 



If 7 = 1, we obtain 
IIE'^+HL 



< 



^ ^io,jo,k+iyio,io io,jo,k+iyjo,joJ io,. 



JO 



/3 fc+1 



«0 JO 



< I IE*' 1 1 no. 



m.=0,m.j^j 



□ 



c?.^. Stability for two-dimensional explicit difference scheme 



Theorem 3.6. //O < 7^ 



< 



1-2-T 



where 



(Ai)" ' (Ay)/3 - r(2-7)C, 



max 

0<i<Af^,0<i<7Vy,0<fc<7Vt 



r(4 - a) 



r(4-/3) 



i/ien i/ie explicit difference scheme \2.19\) of the two-dimensional time-space Caputo-Riesz 
fractional diffusion equation U.l\) with 0<7<1, 1 < a, (3 < 2 is stable. 

Proof. Let • (z = 0, 1, ... , A^^; j = 0, 1, ... , Ny-, k = 0, 1, ... , Nf) be the approximate 



solution of u^ j, which is the exact solution of the explicit scheme fl2.19p . Denoting e 



14 



'fj — u^j and using f l2.19p . we obtain the following perturbation form 



~ (1 + '^i,j,09i,i + '^i,jfi9j,j)^i,j + '^i,j,0 9i,m^m,j + '^i,j,0 9j,m^i,my — 0, 

m=0,mj^i m=0,m^j 
k-1 

(1 - /i + al^^,g^, + al,,k9iM,j + J^^h - h^MT + 



When 7 = 1, (13.41) can be rewritten as 



(3.4) 



Ny 



Using the mathematical induction, we can prove the desired result. Under the conditions 

k IT 



of the theorem, there exists 1 + a'. , .g'^ , + > 0. Let u'^ = [u^, u^, . . . , u 



where vq = [efo) ^i,!; • • • > ^i^ATy] and ||E ||oo = max IhjI- For k = 0, supposing 



lei . I 



I El 



max \e}A, we obtain 



0<i<N^,0<j<Ny 



0<i<N^fl<j<Ny 



"^O JO I 



Ny 



\(1 _|_ n-' -I- rr" V° -I- rr' p° -\- rr" \^ p° I 

m=0,mj^io m=0,mj^jo 



K (1 -\- rr' -\- rr" )\f^ \ -\- rr' \ ^ i o i i « \ ^ ^ i o 

m=0,m^io m=0,mj^jo 



N^ 



Ny 



IE" 



rr' ST' n°' llpOll -L rr" \^n'^ 

"io,jofi/_^9io,m\\^ lloo Ojgjg^o^ yjo.mll-^ 



< (1 + <,Jo,0<,^o + K,,j0,09, 
= ||E°||oo + O-io,jo,0 X^5'i'o,ml|E°||oo + CTiljo.O X] 5'jo,m 1 1 E° 1 1 oo < ||E' 



m=0 



m=0 



Assuming ||E^|U < \ \E'\U k = 1,2, . . . , k, and |ej+i| = ||E^-+i|| 



we get 
I lE^'+^l 



max 

0<i<N:c,0<j<Ny 



I fc+l| 



I fc+1 I 



fe-1 



s=l 



N^ 



Ny 



I / a k I II /3 fc 

io,jo,k / J 9iQ.m^m,jo ' ^io,jo,k / j 9 jQ.m^io,m 



15 



fc-1 

< (1 - /l + (ri^^j^^k9Z,io + '^'Ljo,k9jo,jo'>\^io,jo\ + ^ ^''+1)14 Jo I + ^^I^LjoI 



/ X ^ a \ k \ II ^ ^ /3 I 

m=0,m^io m=0,mj^jo 
fc-1 



fc I 
io,m I 



< (1 - /i)||E'=|U + J2ils - Wi)||E'=-^|U + /fc||E°|U 

s=l 

Ny 

-4- rr' \^ IIF'^II -4- rr" \^ IIF'^ 

"io,jo,k 2^ yicmll-^ Moo -r ^io,jo,k yjcmll-^ 

m=0 m=0 

0| 



< (l-/i)||E°|U + 5^(/,-Wi)||E°|U + /fc||E' 

s=l 

fc-1 

= (l-/i + $^(/.-Wi) + /fc)l|E°||oo = |IE^ 
=1 



00 



When 7 = 1, it is easy to check that 



i-pfc+ii I _ I ,fc+i I 
1-^ Moo Fiojol 



fcl 



1 00 



fc 

— i'^ ~^ '^io,jo,k9io,io + '^io Jo.fc^iojo) 1^*0 JO I ~'~ ^io,jo,k 5'io,mkm,jol ^k)jo,k 9jo,m\^io,' 

m=0,mjtio m=0,m^jo 

Ny 

^ (1 + ^io,jo,k9io,io + '^iojo,fc5'fo,jo)ll^''ll°o+ ^«'o,io,fc 5'i°,ml|E''||oo+ "^iojo.fc 5'jo,mll^ 

Ny 

= I IE''! loo + '^io,jo,k X/ 5'i^,ml|E''|U + <^io,jo,k 5'jo,ml|E''||oo < HE^'Hoo < ||E°||oo- 

m=0 m=0 

□ 

4. Convergence analysis 

We use four subsections to prove that the global truncation error of the schemes f l2.13p 
and (12.181) used to solve f l2.10p is 0{t'^~'^ + (Ax)^), and the global truncation error of the 
schemes (|2A6|) and fl239|) used to solve ([H]) is C(r2"^ + {Axf + (Ay)^). 

4.1. Convergence for one- dimensional implicit difference scheme 

Theorem 4.1. Let he the approximation solution ofu{xi,tk) computed by use of the 
implicit difference scheme 112. 13\) . then there is a positive constant C such that 

\uixi, tk) -u^\< C{r^-^ + (AxY), i = 0,l,...,N,;k = 0,l,...,Nt. 
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Proof. Let u{xi,tk) be the exact solution of fl2.10p at the mesh point {xi,tk). Define 
= u{xi,tk) — u^, and = [eQ,e\, . . . Subtracting (12.111) from (12.121) and using 

e° = 0, we obtain 



m=0,m^i 

k-1 



m=0,m^i s=0 

When 7 = 1, (14.11) can be written as 

(1 - ^i^&l - Stm^^m = Rl k = 0, 

m=0,mj^i 

(1 - u,Mi&'^' - ^^Ml E ' = + ^ > 0- 



(4.1) 



(4.2) 



Denoting that 1 16*^1100 = max \e^\ and -Rmax = max \Rf\, the desired result can 

0<i<Nx 0<i<N^,0<k<Nt 

be proved by using mathematical induction. 

(1) Case < 7 < 1: For k = 0, supposing \e}\ = ||e"^||oo = max \e}\, according to 

" 0<i<N^ 

Lemma 3.2, we get 



m=0 m=0,m=iiQ 
I RiQ I — Riaax ^0 Riaax • 



Supposing ||e''||oo < /r\-Rmax, k = l,2,...,k, and jef+^l = ||e^+^||oo = max k- 



_ ii^fc+i|| _ i^fc+i| 

according to the result of Lemma 3.1, < k = 0,1,..., A;, therefore, He'^Hoo < 
l^^Rmsux, k = 1,2, . . . , k. Then we have 

A'. 

a |^fc+l| 



F lloo - l^io I S l^io I -Wio,A:+l / , 9io,m\^i. 



«0 

m=0 

A;+l| 



Af, 

< (1 - '^io,fc+l5'i'o,io)kir^l -^i(),fc+l E ^'i'o.mkm'^^l 



*0 
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< 



(1 - ^io,fc+lC,io)4^^ 
k—s 



k+1 



k+l 



k-1 



fc+1 



s=0 



k-l 



\{i-h)ei + Y.{is-is^M;' + K 



k+l 

io 



s=l 



From Lemma 3.1, we obtain 



fc-i 



\^ oo 



< (l-/i)||e'=||oo + $^(/.-/ 

s=l 

/ k-l 

<ik'\^-h + Y.iis-i 



:+i)l|e'"1loo 



\R 



io I 



(2) Case 7 = 1: Using similar idea leads to 



,fc+ii 



m I \^io' 



< I|e1|oo + K^'| < (fc + l)i?n.a. < C(A; + l)r(r + (Ax)2) 

< C"(r + (Ax)') = C\t^-^ + (Ax)2). 



fc+i| 



□ 



Convergence for one- dimensional explicit difference scheme 
Theorem 4.2. Let he the approximation solution ofu{xi,tk) computed by use of the 



< 



r(4-a)(l-2-T) 



explicit difference scheme ^KW- U < ^ - A.^Cma.r{2-^){i-2^-) 

max cixiAu), then there is a positive constant C such that 

0<i<Na:,0<k<Nt 



where C„ 



\uixi, tk) -u\\< C(r2-^ + (Aa;)2), z = 0, 1, . . . , iV,; A; = 0, 1 



; • • • ; 



Proof. Define = u{xi,tk) — and i?max = max \Ri\- Analogously, we have 

0<i<Nx,0<k<Nt 



r] , 



k = 0. 



k-i ('43) 
= (1 - /i + a,^kgti)e1 + J2^ls - h+i)et' + ^t,k atm^t + R^^\ k > 0. 



s=l 
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When 7 = 1, (14.31) can be rewrote as 
e] = Rl k = 0, 

^ (4.4) 

ef+i = (1 + a.,kgti)e1 + a.,k J2 dtmsl + R^^\ k > 0. 

m=0,mj^i 

Denoting that = \e'^,e'^, . . . ,e\,] and ||e'^||oo = max lef^l, we use mathematical induc- 
tion to prove the desired result. 

(1) Case < 7 < 1: For k = 0, supposing \e}, I = ||e^||oo = max \e}\, we have 

" 0<i<N^ 
11^ lloo I^JqI I'^ifll — "^max ^0 -^max- 

Assuming He'^Hoo < /f\-Rmax, k = 1,2, . . . ,k, and |e(;+^| = He'^+^Hoo = max |ef+^|, using 
^~ ^ < ^k'^, k = 0,1, . . . , k, therefore, | le'^l |oo < lk^Rma.x, k = 1,2, ... ,k. Then we get 



|pfc+l|| _ ufc+i 
F Moo 



s=l m=0,mj^io 
k-1 

< (1 - /i + ^^o,k9Z,oM\ + J2(^s - + ^^o,k E + 

s=l m=0,m^io 

< (1 - /Olle'^IU + - Wl)l|e^"1loo + ^io.fc Xl^*"o,ml|e1|oo + i?max 

s=l m=0 



< ^fc ""^ 1^1 — ^1 + ^^(^s ~ ^s+l) + 4j -Rmax 
= /,T^i?max<C(r2"^ + (Ax)2). 

(2) Case 7 = 1: Analogously, we have 

lle'+'lloo = I4^'| < (1 + a^oAMl + a.„k 9lJei\ + R„ 

— 1 1^ I loo ~l" ^«o,fc ^ ^ 5'joi™-l 1^ I loo ~l" Rn 



"rnax 



m=0 

< ||e1|oo + i?max<C(r2-^ + (Aa;)2). 

□ 

^.c?. Convergence for two-dimensional implicit difference scheme 

Theorem 4.3. Letu^j be the approximation solution of u{xi,yj,tk) computed by use of 
the implicit difference scheme ^2.10^ . then there is a positive constant C such that 

\u{xi,yj,t,) - uU < Cir'-^ + (Ax)' + (Ay)'), 
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where i = 0,1, . . . , N^; j = 0,1, . . . , Ny] k = 0,1, . . . , Nt. 

Proof. Defining e^j = u{xi,yj,tk) 
ing fl2.14p from fl2.15p . we obtain 



Proof. Defining , = u{xi,yj,tu) — uf,, Rmax = max \Rf \, and subtract- 

V ^/ Jj' max o<j<7V,,0<i<Afj;,0<fc<Aft ' ^-^ ' 



(1 ~ ^i,j,k+l9i',i ~ ^i,j,k+l9jj)^i,j ~ ^i,j,k+l 9i',m^m,j ~ ^i,j,k+l 9j,m^i,m 

= R\j, k = 0, 

l-"- ^i,j,k+l9i,i ^i,j,k+l9jj)^i,j ^i,j,k+l 2_^9i,m'^m,j ^i.i,k+l / j 9i,m'^i,m 

= 5^(/.-WiKT^ + i?f;f, k>o. 

s=0 

When 7 = 1, (14.51) can be written as 



(4.5) 



A^. Ny 
l-^ ^i,j,k+iyi,i ^i,j,k+iyj,j)^i,j ^i.j.fc+l / ^ ■tft.m^m.j '^ij.fc+l / ^ iJj.m'^i.ni 

= el^ + R^+\ k>0. 

Denoting u^' = [ug, , . . . , u^J^, uf = [ef q, i, . . . , ef^^J^ and | |e*^| U = ^^.^ max \e 
we prove the desired result by mathematical induction. 



(1) Case < 7 < 1: For k = 0, supposing \e}^,\ = ||e^||oo = max kKL we 

^ ' 5 r-r- o I lojQ 0<i<N^,0<j<Ny' 



have 



e 



1 



^1 I < 1^1 \ -m' \^ Ir^ I - Ld" 'S^ If ^ 



«o JO I 



m=0 m=0 

< 



^io,joA9io,io ^io,joA9jo,j,i)^io,jo 

Ny 

' I "^iQ Jo I — -^max ^0 -fi/max • 



fc+ll 
« J I ' 



Assuming ||e''||oo < ^r^Rmax, k = 1,2, . . . ,k, and le^^ll = He'^+^Hoo = max |e- 

using < l'^^ , k = 0,1, ... ,k, therefore, ||e^'||oo < ""^Rmax, k = 1,2, ... ,k. Then we 
have 



|pfc+l|| = Ufc+l I < ufc+l I _ / \f:>'+^\—i," \^ n'^ If'' 



iojo I 



m=0 m=0 



«o,jo,fc+iyjoj'ol «o,iol '^Jojo.fc+l / ^ "io.ml^iojol 
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<C (1 — rj' — uj" Mf 



E a I fc+1 I 

m=0,m^io 



m=0,m^jro 



< 



(]_,,' —I )" rr \f^^^ 



— u: 



io JO, ^ ^ 3- 

m=0,m^io 



a fc+1 _ // 
io,m^m,jo ^io,io,fc+l 



fc-l 



m=0,m7^jo 



/3 fc+1 
jo,m io,m 



< 



< 



s=0 



^k—s I 



fc-1 



^fc ' 1 - ^1 + J2(^s - 1 



+ly 



(2) Case 7 = 1: Using similar idea leads to 



■>fc+i| 



< 



V '^iojo,fc+l"«o,jo «o,jo,fc+l"ioj'o'' «oj'o 

AT, iV. 



^fc+1 



y 



«o Jo, fc+1 / > fc+1 



i4.o + ^Sl<ll^1loo + l^Sl<(^ + i)R. 



<C(r^-^ + (Ax)^ + (Ay)^). 



4.4- Convergence for two-dimensional explicit difference scheme 

Theorem 4.4. Let u^j be the approximation solution of u{xi,yj,tk) computed by use 
the implicit difference scheme l\2.19\) . If < + ^ r(2-'f)C — ' where 

o<i<7v,,o</<7v„o<fc<7Vt I r(4-a) ' r(4-/3) 
t/ien t/iere zs a positive constant C such that 

\u{xuyj,t,) - ul^\ < C{r'-^ + (Ax)2 + (Ay)'), 
where i = 0,1, . . . , N^; j = 0,1, . . . , Ny] k = 0,1, . . . , Nt. 
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Proof. Define ef, = uixi^VjAu) —uf, and Rmax = max \RfA- Similarly, 

'■^ V J / ^,3 0<i<N^,0<j<Ny,0<k<Nt 

we obtain the following form 



k-l 



s=l 

Ny 



(4.6) 



m=0,mj^i m=0,m^j 



When 7 = 1, (14.61) becomes 

^i.i — ^i.i^ k = 0, 



^ij^ "(-'- + ^i,j,k9?,i + '^i,j,k9j,j)^i,j 

Ny 

Using mathematical induction, we can obtain the desired result. Letu'^ = [uq, u^, . . . , u^^]-^, 
where = [£^(,,£^1, . . . »f 1^ and lle^'lloo = max \efA. 

(1) Case < 7 < 1: For k = 0, supposing \£},j„\ = ||e^||oo = max ki^L we 



obtain 



fc+l| 



Assuming He'^Hoo < /r \Rmax, k = 1,2, . . . ,k, and jef+M = He^'+^Hoo = max |e- 

using /-""^ < Z^"*^, = 0,1, ... ,k, therefore, He'^Hoo < ""^Rmax, k = 1,2, ... ,k. Then we 
have 



fc+l|| _ \^k+l 



|G I loo 



'*0 JO I 



k-l 



- (1 /i + cr[^j^^k9io, to + (^io,jo,k9jo,jo^^io,jo + X]*^^^ 

iV:. Ny 

I ^' „a k I // B k I p, 

k-l 

< (1 - /l + (^'io,jo,k9?o,io + ^'Ljo,k9jo,jo'>\^ioJo\ + X]'^^'* ~ ^'^+^)\^io,jo 

s=l 

Ny 

+ ^iojo.fc fi'io,mkmjol ^«OJO,fc 5'jo,mkio,ml + R-max 

k-l 

<(l-/i)||e'=||oo + 5^(/.-Wi)||e'=-^' 
=1 
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'^io,jo,k'^^ 9io,m\\^ l|oo ^in.in.fc ^^ ^fn-mll^ Hoo + R-max 
m=0 m=0 

-1 ^ 



= /fc 'Rrnax < C (t''-^ + (Ax)^ + (Ay)^). 

(2) Case 7 = 1: It is easy to check that 



1*^ Moo \^io,jo\ 



max 



— ^io,jo,k9io,io ~^ ^io,jo,k9jo,jo)\^io,jo\ 
_|_ ^/ „a \ k I I // B I k I _|_ -R 

— + ^io,jo,k9io,io + '^iljo,fcfi'io,io)ll*^''ll°o 

+ ^io,jo,k ^ ] 5'io,ml 1^ I |oo+ (^iojo^k ^ ] 5'jo,ml 1^ I |oo + R-max 

— 1 1 ^ 1 1 oo ~l~ Rmax — i^k -\- 1 ) Rjnax 

<C(r2-^ + (Ax)2 + (Ay)2). 

□ 

Remark 4.4. Let u G (0,1), there exist xlDxU = DxJ^x^^~''^u = D'^^l^x^'^'^^'u, and 
xD^^u = —DxDxji^ '^^u = D'^xDx^ ''^u [sl- So after carefully dealing with the boundary 
conditions, all the above presented numerical schemes still work well when the order of 
space fractional derivatives a G (0,1) and/or /3 G (0,1). And the convergence rates 
remain, and all the theoretical analyses are valid. Some of the numerical results are also 
given in Table HI 



5. Numerical results 

In this section, we numerically verify the above theoretical results including conver- 
gence rates and numerical stability. And the /qo norm is used to measure the numerical 
errors. 

5.1. Numerical results for the implicit scheme for the one- dimensional time-space Caputo- 
Riesz fractional diffusion equation 

Consider the one-dimensional time-space Caputo-Riesz fractional diffusion equation 
fl2.10p . on a finite domain 0<x<l, 0<t< 1/2, with the coefficient c{x,t) = x°'t^~'' , 
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the forcing function 



f{x,t) =-r(3 + 7)tV(x-l)^ 



3™Q 



+ 



cos{a7r/2) 



X 



2-a 



:i-x 



,2-a 



-6- 



X 



3-a 



+ (l-x 



|3-Q 



12- 



X 



4— a 



+ (l-x 



A— a 



r(3-a) " r(4-a) 1(5 - a) 

the initial condition u{x, 0) = 0, and the boundary conditions u{0,t) = u{l,t) = 0. This 
fractional diffusion equation has the exact value u{x,t) = t'^^^x'^{l — x)^, which may be 
confirmed by applying the fractional differential equations 



d:(x - x,r 



r(t; + l) 
r(t; + 1 - u) 



[X - XL 



,D^^{xr-x) =— — — -{xb-x) 



r t; + 1 - u) 



Table 1: The maximum errors and convergence rates for the implicit scheme fl2.13p of 
the one-dimensional time-space Caputo-Riesz fractional diffusion equation fl2.10p with 
variable coefficient c{x,t) = x°'t^~'' at t = 1/2, and the time and space stepsizes are 
equal, i.e, r = Ax. 



T, Ax 


a = 1.2,7 = 0.9 


Rate 


a = 1.2,7 = 0.5 


Rate 


a = 1.2,7 = 0.1 


Rate 


1/40 


3.1438e-004 




6.3187e-005 




2.7395e-005 




1/80 


1.4713e-004 


1.0954 


2.2183e-005 


1.5102 


7.6703e-006 


1.8366 


1/160 


6.8748e-005 


1.0977 


7.7888e-006 


1.5100 


2.0317e-006 


1.9166 


1/320 


3.2097e-005 


1.0989 


2.7378e-006 


1.5084 


5.2912e-007 


1.9410 




a = 1.9,7 = 0.9 


Rate 


a = 1.9,7 = 0.5 


Rate 


a = 1.9,7 = 0.1 


Rate 


1/40 


2.7655e-004 




5.6774e-005 




2.1114e-005 




1/80 


1.2919e-004 


1.0981 


1.9669e-005 


1.5293 


5.4717e-006 


1.9482 


1/160 


6.0294e-005 


1.0994 


6.8351e-006 


1.5249 


1.4145e-006 


1.9517 


1/320 


2.8133e-005 


1.0997 


2.3841e-006 


1.5195 


3.6518e-007 


1.9536 



Table [T] shows the maximum errors, at time t = 1/2 with r = Ax, between the exact 
analytical values and the numerical values obtained by applying the implicit scheme f l2.13p . 
Since the scheme has the global truncation error (9(r^~'^ + (Ax)^), the convergent rate 
should be 2 — 7 being confirmed by the numerical results. 

Table |2] shows the maximum errors at time t = 1/2, and the time and space stepsizes 

2 

are taken as r = (Ax) 2-7 . The numerical results confirm the 2nd order convergence 
in space directions. In particular, the numerical results when a = 0.3 G (0, 1) are also 
presented, which confirm the statement of Remark 14.41 
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Table 2: The maximum errors and convergent rates for the imphcit scheme f l2.13p of the 
one-dimensional time-space Caputo-Riesz fractional diffusion equation (12.101) at t = 1/2 

2 

with variable coefficient c(x, t) = x'^t and r = (Ax)^-^, where 7 = 0.9. 



T, Ax 


a = 1.9 


Rate 


Q = 1.5 


Rate 


a = 1.2 


Rate 


a = 0.3 


Rate 


1/10 


2.4699e-004 




2.5801e-004 




2.5510e-004 




2.4135e-004 




1/20 


6.2966e-005 


1.9718 


6.5569e-005 


1.9763 


6.4560e-005 


1.9823 


6.1043e-005 


1.9832 


1/40 


1.5697e-005 


2.0041 


1.6226e-005 


2.0148 


1.5934e-005 


2.0185 


1.5085e-005 


2.0167 


1/80 


3.9368e-006 


1.9954 


4.0475e-006 


2.0032 


4.0433e-006 


1.9785 


3.7583e-006 


2.0050 



5.2. Implicit schemes results for two-dimensional time-space Caputo-Riesz fractional dif- 
fusion equation 

We further examine the two-dimensional time-space Caputo-Riesz fractional diffusion 
equation (11 .ip . on a finite domain 0<a;<l, 0<t< 1/2, with the variable coefficients 
c{x,y,t) = 2x'^y^t^~"' , d{x,y,t) = 2x^y'^t^~"', the forcing function 

f{x,y,t) 

= lr(3 + 7)tV(x-i)V(y-i)' 



2 



+ 
+ 



cos{a7i/2) 

t^X^+f^ix - 1)2?/" 

COS (/37r/2) 



x^-" + (1 - ^x^-" + (1 + (1 -x)^-" 
6 —, ^ h 12- 



r(3-a) r(4-a) r(5 - a) 

y'-' + (1 - y?-' , y'-^ + {i-y?~' ^ ,^y'~' + (i - y)'-^ 

O z— — h 12 



r(3-/3) r(4-/3) r(5-/3) 

the initial condition UQ{x,y) = 0, and the boundary conditions B{x,y,t) = 0. It has the 
exact solution 

u{x,y,t) = t'+^x'{l-x)Y{l-yy- 

Table E] shows the maximum errors of the implicit scheme (I2.16P at time t = 1/2, 
and r = Ax = Ay. Since the implicit scheme (I2.16P has the global truncation error 
C(r2~'^ + (Ax)2 + (AyY), the convergence rate should be 2 — 7 being confirmed by the 
numerical results. 



6. Conclusions 

This paper discusses the finite difference schemes for the time-space Caputo-Riesz 
fractional PDEs with variable coefficients, and the order of time fractional derivative 
belongs to (0, 1) and the order of space fractional derivatives locate in (1, 2) . The obtained 
schemes have (2 — 7)-th order convergence rate in time and 2nd order convergent rate 



25 



Table 3: The maximum errors and convergence rates for the imphcit scheme fl2.16p of the 
two-dimensional time-space Caputo-Riesz fractional diffusion equation (11 .ip at t = 1/2 

and r = Aa; = Ay, with variable coefficients c(x, y, t) = 2x°'y'^t^~"' and d{x, y, t) = 
2x^yn^-^. 



T, Ax, Ay 


7 = 0.9 




7 = 0.5 




7 = 0.1 






a = 1.2,^ = 1.3 


Rate 


a = 1.2,^ = 1.3 


Rate 


Q = 1.2,13 = 1.3 


Rate 


1/10 


7.7867E-005 




3.1936E-005 




2.2077E-005 




1/20 


3.6381E-005 


1.0978 


1.0357E-005 


1.6246 


5.6507E-006 


1.9660 


1/30 


2.3084E-005 


1.1220 


5.4951E-006 


1.5632 


2.5497E-006 


1.9627 


1/40 


1.6839E-005 


1.0965 


3.4925E-006 


1.5755 


1.4509E-006 


1.9598 




a = 1.8,^ = 1.7 


Rate 


a = 1.8,13 = 1.7 


Rate 


Q = 1.8,13 = 1.7 


Rate 


1/10 


7.8209E-005 




3.1251E-005 




2.0743E-005 




1/20 


3.6912E-005 


1.0832 


1.0369E-005 


1.5917 


5.5573E-006 


1.9002 


1/30 


2.3535E-005 


1.1099 


5.5386E-006 


1.5465 


2.5115E-006 


1.9588 


1/40 


1.7122E-005 


1.1060 


3.5433E-006 


1.5527 


1.4310E-006 


1.9552 



in space. The detailed numerical stability analysis and error estimates are presented, 
and the extensive numerical experiments are performed, which confirm the theoretical 
results. In particular, the numerical schemes still work well for the time-space Caputo- 
Riesz fractional PDEs with the order of its space derivatives belongs to (0, 1), and all the 
theoretical analyses are still valid. 
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